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Background 

•  A  vector  sensor  employs  multiple  (Nc)  sensor  components.  For  acoustic: 
Nc  =  4;  for  electromagnetic:  Nc  =  6. 

•  Two  strategies:  matrix-based  scheme  (left)  and  tensor-based  scheme  (right). 


•  Tensor-based  scheme  conserves  the  multidimensional  structures  of  the  data, 
and  provides  better  performance  than  the  matrix-based  scheme. 

•  Our  goal:  Apply  the  nested-array  strategy  to  vector-sensor  arrays  via  tensor 
modeling. 

•  Challenge:  Multidimensional  operation. 


N  X  Nr 


NNC  x  1 


■  III 


I 


5 


NNC  x  T 


Concatenation 
of  components 

Multiple 

samples 


N:  Sensors 
T :  Snapshots 
Nc:  Components 
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Signal  Model 

•  Consider  a  two-level  nested  linear  array  with  N  vector  sensors, 
each  vector  sensor  is  an  iVc-dimensional  vector. 

•  Assume  K  far-field  sources  from  directions  {(</>&,  #&),  k  =  1, 
q f>k  and  9k  represent  the  azimuth  and  elevation  angles 

•  Then,  we  obtain  the  tensor  measurement  model: 


Nc 

Y ( t )  =  *A,  X3  x(t)  +  E(t ), 


•  Y(t):  N  x  Nc  measurement  matrix  at  time  t 

•  A:  N  x  Nc  x  K  array  manifold  tensor 

•  x(t):  K  x  1  source  signals 

•  E(t):  N  x  Nc  measurement  noise  at  time  t 

Ax3B:  Mode-3  product  of  A  G  C/lx72x73  and  Be  CJ1X  J2x/3,  defined  as  (Ax3  = 


The  output  of 
. . . ,  K},  where 


(i) 

^3  aH*2*3&ili2*3 
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Signal  Model  (Cont.) 

,j,k  —  where  Aj~  —  &kPk 

dk  =  [ei27TUkri/x, . . . ,  e^27TUkrN/x]T  is  the  phase  delay  vector. 

Uk  =  [cos^cos^/e,  sin0fccos0fc,  sm0k]T  is  the  unit  vector  at  the  sensor  pointing 
towards  the  fcth  signal. 

Pk  is  the  steering  vector  of  a  single  vector  sensor  located  at  the  origin. 

Based  on  (1),  we  get  the  N  x  Nc  x  N  x  Nc  interspectral  tensor: 

n  =  E [Y  O  Y*]  =  A  x 3  Rxx SA*  +  E [E  o  £7*],  (2) 


which  is  a  tensor  version  of  the  covariance  matrix  in  the  scalar  case. 
We  apply  mode-2  matricization  to  7 Z,  and  obtain 

Z  4  Rf  =  {A*T  @A)x3S  +  aff  5 


-  Z  is  an  NCN2  x  iVc  matrix 

_  s  _  r  2  2  IT 

*  —  [cr1 ,  a2,  .  .  .  , 


-  =  blkdiag(Tf . . . ,  if),  where  =  [ef ,  e^, . . . ,  e^]T 


(3) 
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Signal  Model  (Cont.) 

Generalized  Kahtri-Rao  product  © 


j 


Consider  one  column  vector: 

04^  ®  Jl\hk  =  (A^\k  O  cA.  jk 

N2Nc  x  1  NNC  x  1  N  xl 

0^(3)  © 

N2Nc  x  Nc  x  K 


*Ax323:  Mode-3  inner  product  of  A  £  C^l  x^2  x^3  and  23  £  C^l  x^2x^3t  defined  as 

A  ©  23:  Generalized  Khatri-Rao  product  of  A  £  C^x^3  and  23  £  C^1x^2x-^3i  with  d  imension  I-^J  x  I2 
(A  ©  e)(f1  +  (j-1)/1),i2,i3  =  aj,i3bii,i2,is 

Ao  B.  Outer  product  of  A  £  C^lx^2  and  B  £  C^1x^2t  defined  as  (A  o  =  a*l*2^1^2 

^2):  mode-2  matrix  unfolding  of  tensor  A  £  C^l  x  J2  x  J3 ,  defined  as  (•A(2))i2,(ii  — l)/3+i3  =  (*^-)f1f2f3 


°'ili2i3bhj2i3 
x  /3,  defined  as 
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Signal  Model  (Cont.) 


•  Comparing  (3)  with  the  original  signal  model  (1),  we  observe  that  Z  in  (3) 
behaves  like  a  received  measurement  with  a  longer  vector-sensor  array  whose 
manifold  is  given  by  © A In  (3)  the  equivalent  source  signal  vector  is 

represented  by  s  and  the  noise  becomes  a  deterministic  matrix  given  by  aft . 

•  Looking  at  the  structure  of  the  tensor  A*^  ©  A,  we  observe  that  there  are 
Nc  sets  of  horizontal  slices,  and  they  have  same  horizontal  slices  except  for 
different  amplitudes.  Without  loss  of  generality  we  will  consider  to  use  only 
the  first  set. 

•  Note  that,  although  we  choose  only  one  set,  we  use  the  information  from  all 
the  Nc  sensor  components. 
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Signal  Model  (Cont.) 


Dimension  change  of  the  manifold  tensor  using  the  nested-array  strategy. 


Manifold  tensor  *4  of  the 
original  signal  model. 
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Signal  Model  (Cont.) 


Dimension  change  of  the  manifold  tensor  using  the  nested-array  strategy. 


After  matricizing  the 
interspectral  tensor  Jl 


•  Nc  parallel  sets  of 
horizontal  slices 

•  Different  amplitudes 
but  same  information 


CSSIP  Lab 


8 


0  Washington 

University  in  StLouis 


Signal  Model  (Cont.) 

Dimension  change  of  the  manifold  tensor  using  the  nested-array  strategy. 


NCN2 


A&  © 


Concentrate 
on  one  set 


Multiple  repeated  horizontal 
slices  (matrices) 
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Signal  Model  (Cont.) 


Dimension  change  of  the  manifold  tensor  using  the  nested-array  strategy. 


Remove  repeated  horizontal 
slices  and  resort  them 
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Spatial  Smoothing 

•  Spatial  smoothing  is  used  to  exploit  the  increased  DOFs  by  building  up  the 
rank  of  the  observation  tensor. 

•  Similar  to  the  scalar  sensor  case,  we  remove  repeated  slices  and  sort  them 
according  to  virtual  sensor  positions: 

Z  —  A,  X  3  «S  ~b  (J2E. 

•  Divide  these  2 N  —  1  (N  =  N2/ 4  +  N/2)  virtual  sensors  into  N  partially 
overlapping  subarrays,  and  the  Zth  subarray  is  Zi  =  Ai  xss  +  XeEl. 

•  Define  TZi  =  Z\  o  Z\  ,  and  take  the  average  of  7 Zi  over  all  l:  ^  YliLi 

•  T  is  the  N  x  Nc  x  N  x  Nc  spatially  smoothed  interspectral  tensor. 

•  Based  on  7~,  the  nested  array  with  N  vector  sensors  provides  N  —  1  DOFs, 
whereas  a  ULA  with  the  same  number  of  vector  sensors  can  provide  only  N  —  1 
DOFs. 
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Higher-Order  Singular  Value  Decomposition  (HOSVD) 

•  HOSVD  is  a  higher-order  generalization  of  the  matrix  singular  value 
decomposition  (SVD)  [3]. 

•  The  HOSVD  of  tensor  'T  can  be  written  as 

T  =  Kxi  U\  X2  U2  X3  Us  X4  L^4,  (4) 

where  Ui,Us  G  C^xi\  and  (72,  U4  G  C NcXNc  are  orthonormal  matrices, 
provided  by  the  SVD  of  the  i-mode  matricization  of  the  tensor  77  7 ~(i)  = 
UiAiV/1 .  K,  G  C^xivcxivxivc  js  ^ e  core  ^-ensor 

•  Since  T  is  an  Hermitian  tensor,  i.e.,  th^i3M  =  h  i 2,  23^4,  the 

HOSVD  of  T  can  be  written  as 

T  =  Kx1U1  x2U2  x3(7*  x4[/2*.  (5) 


[3]  M.  Boizard,  G.  Ginolhac,  F.  Pascal,  S.  Miron,  and  P.  Forster,  Numerical  performance  of  a  tensor  MUSIC  algorithm  based  on  HOSVD  for 
a  mixture  of  polarized  sources,  in  EUSIPCO  2013,  Marrakech,  Marocco,  Sep.  2013. 
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Applications:  Acoustic  Vector  Sensors 

•  For  acoustic  vector  sensors,  Nc  =  4. 

•  The  array  steering  matrix  can  be  written  as  Ak  =  dkP^i  with 


Pk  =  [1  ,U 


TiT 
k\  5 


which  is  the  steering  vector  of  a  single  vector  sensor  located  at  the  origin. 
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Applications:  Electromagnetic  Vector  Sensors 

For  EM  vector  sensors,  Nc  =  6.  Here,  we  consider  polarized  signals. 

The  array  steering  matrix  can  be  written  as  Ak  =  dkP with  pk  =  Vkpk- 
where 


V 


— sin0fe 

— cos^sin  6k 

COS  (j)k 

— sin^fesin^fe 

0 

cos  9k 

cos^fesin^fc 

sin  4>k 

■sin(pksindk 

-cos  <t>k 

cos  9k 

0 

pk  =  [cos7fc  sin7fee-?%]T. 

and 


J 


Ak  is  the  N  x  Nc  steering  matrix  of  the  array  associated  with  a  polarized 
signal  coming  from  the  direction  (</>&,  #&)  with  polarization  (7 k,Vk),  where 
7 k  E  [0,27 r]  and  rjk  E  (— 7r,  tt]  are  polarization  parameters.  T4  is  the  steering 
matrix  of  one  EM  vector  sensor  associated  with  the  fcth  signal,  pk  is  the 
polarization  vector  for  the  kth  signal. 
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Acoustic  Case  I:  MUSIC  Spectrum 


9  (rad) 


Fig.  1:  MUSIC  spectrum  using  a  nested  acoustic  vector-sensor  array  with  6  sensors,  as  a  function 
of  elevation  angle  0,  T  =  1000,  SNR  =  OdB,  K  =  6  sources  with  the  same  azimuth  angles. 

Observation:  The  2-level  nested  array  resolves  all  the  6  sources. 
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EM  Case  II:  DOA  Estimation 


Fig.  2:  MUSIC  spectrum  using  a  ULA  (left:  tensor- based)  and  a  nested  array  (middle:  matrix- 
based;  right:  tensor-based)  with  6  EM  vector  sensors,  as  a  function  of  azimuth  0  and  elevation 
angles  6,T  =  1000,  SNR  =  21.97dB,  K  =  2. 

Observation:  The  proposed  nested  vector-sensor  array  strategy  outperforms  the 
ULA  with  the  same  number  of  sensors.  In  addition,  the  tensor-based  method 
outperforms  the  matrix-based  method. 
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EM  Case  III:  Source  Number  Detection 


Fig.  3:  Probability  of  detection  versus  SNR  using  a  nested  array  with  6  EM  vector  sensors  and 
ULAs  with  6  and  12  EM  vector  sensors,  T  =  1000,  K  =  2. 

Observation:  We  can  see  that  the  detection  performance  of  all  the  three  arrays 
improves  with  increasing  SNRs.  In  addition,  the  nested  array  outperforms  the 
corresponding  ULA  with  same  number  of  sensors  and  performs  closer  to  the  ULA 
with  double  number  of  sensors. 
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Summary 


•  Nested  vector-sensor  array  processing  via  tensor  modeling 

-  We  established  the  analytical  foundation  of  the  nested  vector-sensor  array 
by  exploiting  multilinear  algebra 

-  We  constructed  corresponding  signal  processing  strategies,  and  verified  their 
effectiveness  through  numerical  examples. 
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